Pairing correlations and odd-even staggering in reaction cross sections 

of weakly bound nuclei 
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We investigate the odd-even staggering (OES) in reaction cross sections of weakly bound nuclei 
with a Glauber theory, taking into account the pairing correlation with the Hartree-Fock-Bogoliubov 
(HFB) method. We first discuss the pairing gap in extremely weakly bound nuclei and show that 
the pairing gap persists in the limit of zero separation energy limit even for single-particle orbits 
with the orbital angular momenta / = and 1 = 1. We then introduce the OES parameter defined 
as the second derivative of reaction cross sections with respect to the mass number, and clarify the 
relation between the magnitude of OES and the neutron separation energy. We find that the OES 
parameter increases considerably in the zero separation energy limit for I = and I = f single- 
particle states, while no increase is found for higher angular momentum orbits with e.g., I = 3. We 
point out that the increase of OES parameter is also seen in the experimental reaction cross sections 
for Ne isotopes, which is well accounted for by our calculations. 
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I. INTRODUCTION 

Reaction cross sections an of unstable nuclei provide 
a powerful tool to study the structure of unstable nu- 
clei such as density distribution and deformation 
For instance, a largely extended structure, referred to 
as "halo", of unstable nuclei such as ^^Li ^^Be 
and i''d9(^ Q jjg^g i^een found with such measurements. 
The halo structure is one of the characteristic features of 
weakly bound nuclei, and has attracted lots of attention 
(see Ref. 01 for a recent discovery of halo structure in 
■^^Ne nucleus). 

Experimentally large odd-even staggering (OES) phe- 
nomena have been revealed in reaction cross sections of 
unstable nuclei close to the neutron drip line, e.g., in the 

isotopes 14,15,16c H 18,19,20c |l 28,29,30Ng[| ^ 30,31,32^5 

0, and 36,37,38]N^g Ylei. we have argued that 

the pairing correlations play an essential role in these 
OES. That is, the OES in reaction cross sections is in- 
timately related to the so called pairing anti-halo effect 
discussed in Ref. On the other hand, there have 

been contradictory arguments whether the pairing gap 
disappears jjj] or persists when a nucleus reaches 

at the neutron drip line, i.e., the single-particle energy 
of the last occupied orbit approaches the zero energy. 
If the pairing gap disappeared, the OES effect might be 
either quenched or disappeared completely, unless the de- 
formation parameter is significantly different among the 
neighboring nuclei . 

In this paper, we first discuss the pairing correlations 
close to the zero energy by the Hartree-Fock Bogoliubov 
(HFB) method. We carry out HFB calculations for the 
neutron 3si/2 orbit in ^^Cr changing the separation en- 
ergy in a mean field potential, and examine different defi- 
nitions for an effective pairing gap parameter. This prob- 
lem is also related with the superfluidity of neutron gases 
in the outer crust of neutron stars llSlI . The second mo- 



tivation of this work, in addition to giving the details of 
the analysis in Ref. is to propose a formula to mea- 
sure the odd-even staggering in the reaction cross sec- 
tions. Notice that the OES of the isotope shift of stable 
nuclei has been discussed mostly to clarify the deforma- 
tion changes in odd-mass and even-mass nuclei [l6| . The 
present issue of OES in the reaction cross sections has 
a similar aspect to the previous study in Ref. [l6| in 
one sense, but different in another aspect since it aims at 
studying the existence of the pairing correlation in nuclei 
close to the neutron drip line. 

The paper is organized as follows. In Sec. II, we dis- 
cuss the pairing correlation in neutron-rich nuclei using 
the HFB method. In Sec. HI, we apply a Glauber the- 
ory in order to calculate reaction cross sections. We in- 
troduce the OES parameter for reaction cross sections, 
and discuss it in relation with the pairing correlations in 
weakly bound nuclei. We then summarize the paper in 
Sec. IV. 



II. PAIRING GAP AT NEUTRON DRIP LINE 

In the coordinate space representation, the HFB equa- 
tions read fTMl 
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is the mean-field Hamiltonian, m being the nucleon mass. 
V{r) and A(r) are the mean-field and the pairing poten- 
tials, respectively, and Ei is a quasi-particle energy. Here, 
we have assumed that the nucleon-nucleon interaction is 
a zero range force so that these potentials are local. The 
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upper component of the pair wave function Ui{r) is a 
non-locahzed wave function if the quasi-particle energy 
Ei is larger than the Fermi energy |A|, while the lower 
component Vi{r) is always localized. The pair potential 
A(r) in general has a larger surface diffuseness than the 
mean field potential V{r), and goes beyond it due to the 
non-localized property of the upper component of the 
wave function Ui{r), that is, due to the coupling to the 
continuum spectra [18|. 

In the mean field approximation without the pairing 
correlations {i.e., A(r) = 0), the halo structure originates 
from an occupation of a weakly-bound I = or I = 1 orbit 
by the valence nucleons near the threshold [1^ [2l| . The 
asymptotic behavior of a single particle wave function for 
s wave reads 



tpi{r) exp(-Q!ir), 



(3) 



where Ui is defined as = ^/2^m\ei\JW with the Hartree- 
Fock (HF) energy e^. The mean square radius of this 
wave function is then evaluated as 



HF 
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which will diverge in the limit of vanishing separation 
energy, jeil — >• 0. It has been shown that this divergence 
occurs not only for s wave but also for p wave, although 
the dependence on \si\ is now (r^)HF oc l/-y/|e7f for I = 1 

In contrast, in the presence of the pairing correlations 
{i.e., A(r) 7^ 0), the lower component of the HFB wave 
function, which is relevant to the density distribution, 
behaves as 



Vi{r) oc exp(-/3ir), 
where is given by, 
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using the quasi-particle energy Ei. With the canonical 
(r) in the HFB theory, the quasi-particle en- 
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ergy may be approximately given by 18] 
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In the zero binding hmit, e. 
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and A 0, the asymptotic behavior of the wave function 
Vi{r) is therefore determined by the gap parameter as, 



Vi{r) (x exp 
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The radius of the HFB wave function will then be given 
in the limit of small separation energy as 



(r-^)HFB = 



f r^\v,{r)\^dr 1 

— (X 
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2mA 



(can) 



(9) 



If the gap parameter A^'^'^'^-' stays finite in the zero en- 
ergy limit of e^j^'^'^\ the extremely large extension of a 
halo wave function in the HF field will be reduced sub- 
stantially by the pairing correlations and the root-mean- 
square (rms) radius will not diverge. This is referred to as 
the anti-halo effect due to the pairing correlations px|. 
It was shown in Ref. Q that this is the main reason 
for the observed OES in the reactions cross sections of 
several drip line nuclei. 

In order to study the behavior of the pairing gap in 
weakly bound nuclei, we carry out HFB calculations for 
the neutrons in ^^Cr nucleus. To this end, we use a spher- 
ical Woods-Saxon (WS) potential, 



y{r) = Vof{r) - — — I ■ s, 



dr 



with 



fir) 



1 



l + exp[(r-i?o)/a] 



(10) 



(11) 



for the mean- field potential V{r). Following Ref. [23, 
we take Vq = -51 + 30(iV - Z)/A MeV, i?o = 1.27^1/3 
fm. Vis = -0.71 Vq MeV-fm2, and a=0.67 fm. For the 
HFB calculations, we use a density-dependent contact 
pairing interaction of surface type, with which the pairing 
potential is given by 



A(r) 



1 - 



Po 



Pn{r). 



(12) 



Here, p{r) and Pn{r) are the total particle density and 
the neutron pairing density, respectively, given by 



p{r) = J2 I^^WI'' 
Pn{r) = -'^u*{r)v,{r) 



(13) 
(14) 



We again follow Ref. 
Kair = -420 MeV-fm3 



22| and take po=0.16 fm"^ and 
pair — — -li-u ivj.cv-1111^ with the energy cut off of 50 
MeV above the Fermi energy. In order to construct the 
proton density, we use the same mean-field potential as 
in Eq. ([10]), but with Va = -51 - 30(A^ - Z)/A MeV. 
We also add the Coulomb potential for a uniform charge 
with a radius of Rq. We discretize the continuum spec- 
tra with the box boundary condition. We take the box 
size of i?box=60 fm, and include the angular momentum 
up to Z = 12. Notice that we determine the pairing po- 
tential self-consistently in this model according to Eq. 
([T^. although the Woods-Saxon potential is fixed for 
the mean- field part. The Fermi energy is also determined 
self-consistently according to the condition for the aver- 
age particle number conservation, 



(A) = J p{r)dr = Jj2 hir)\^dr. 



(15) 
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FIG. 1: (Color online) The mean square radii and the paring 
gap for '^^Cr nucleus as a function of the single particle energy 
ews for the 3si/2 orbit in a Woods-Saxon potential. The top 
panel shows the mean square radius of the 3si/2 wave function 
obtained with (the solid line) and without (the dashed Une) 
the pairing correlation. The middle panel shows the root- 
means-square (rms) radii for ^*Cr (the dotted line), ^^Cr (the 
dashed line), and ^^Cr (the solid line). These are obtained 
with the Hartree-Fock ('^''Cr and ^^Cr ) and the Hartree-Fock- 
Bogoliubov C^Cr ) methods. The bottom panel shows the 
effective pairing gap for the 3si/2 state. The solid line is 
the pairing gap defined with the canonical basis, that is, the 
expectation value of the pair potential with respect to the 
canonical basis wave function for the 3si/2 state, while the 
dot-dashed line shows the lowest quasi-particle energy for the 
Si/2 states. 



These self-consistencies are particularly important to in- 
crease the pairing gap for an extremely loosely bound 
orbit [T^. 

The top panel of Fig. 1 shows the mean square ra- 
dius of the 3s]^/2 state. In order to study the dependence 
on the binding energy, we vary the depth of the WS po- 
tential Vq for neutron s wave states while we keep the 
original value for the other angular momentum states. 
We also arbitrary change the single-particle energy for 
the 2d5/2 state from -0.38 McV to -0.05 MeV so that 
the 3si/2 state lies below the 2d5/2 state. The dashed 
line is obtained with the single-particle wave function for 



the mean-field Hamiltonian h, while the solid line is ob- 
tained with the wave function for the canonical basis in 
the HFB calculations. The radius of the single-particle 
wave function for the s-wave state increases rapidly as 
the single-particle energy e approaches zero, and eventu- 
ally diverges in the limit of ews — J" 0. In contrast, the 
HFB wave function shows only a moderate increase of 
the radius even in the limit of ews ~^ 0. The middle 
panel of Fig. 1 shows the root-mean-square (rms) radius 
for the whole nucleus, by taking into account the contri- 
bution of the other orbits as well. For comparison, we 
also show the rms radii for ^*Cr and ^^Cr obtained with 
the same mean-field potential but without including the 
pairing correlation. Notice that the rms radius of ^^Cr 
is larger than that of ^^Cr in the range of ews < —0.56 
MeV. This is due to the coupling of single-particle wave 
functions to a larger model space, including continuum, 
induced by the pairing correlations. On the other hand, 
in the limit of ews ^ 0, the rms radius of ^^Cr shows a 
divergent feature while that of ^^Cr is almost constant, 
for the 3si/2 state. The solid line shows the pairing gap 

evaluated with the canonical basis, a|'^™^ for i — 3si/2j 
while the dot-dashed line is the lowest quasi-particle en- 
ergy Ei. Notice that the right hand side of Eq. (O is 
simply a diagonal component of the HFB Hamiltonian in 
the canonical basis, while the left hand side is obtained 
by diagonalizing the HFB matrix [l^. Due to the off- 
diagonal components, Eq. ^ holds only approximately, 
and thus Ei may be smaller than A^'^'^"^ in actual calcula- 
tions. One can see in the figure that the effective pairing 
gaps persist even in the limit of ews ~^ 0, leading to the 
reduction of the radius of HFB wave function as is shown 
in the top panel of Fig. 1. 

In the simplified HFB model of Ref. [ll|, it was 
claimed that the effective paring gap is diminished or 
quenched substantially for low I orbits with I — and 
1 = 1. In this model, the radial dependence of the 
pairing field A(r') is fixed either as a Fermi-type func- 
tion (volume-type) or a derivative of the Fermi function 
(surface-type) with the same surface diffuseness param- 
eter as in the mean-field potential. Furthermore, the 
Fermi energy is set equal to the single-particle energy, 
e, for the mean-field Hamiltonian /i, that is, A = e. The 
effective pairing gap is defined in Ref. [ill be identical 
to the corresponding quasi-particle energy. If the energy 
in the canonical basis were the same as the single-particle 
energy, Eq. ([7]) indeed yields 

E, = A^^K (16) 

Care must be needed, however, since in general g^'^^"-' 
deviates from the single-particle energy, Ci in the HFB. 
Moreover, when the effective gap is plotted as a function 
of e as has been done in Ref. [ll|, setting A = e leads to 
a violation of particle number in this model, whose effect 
may be large in the limit of e — > 0. 

In order to investigate the consistency of the simplified 
model of Ref. [llj , we repeat the same calculations shown 
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FIG. 2: (Color online) Same as Fig. 1, but with the simplified 
HFB model, in which the pair potential is assumed to be 
proportional to the derivative of the Fermi function. Only the 
Fermi energy A is determined self-consistently in this model. 
The dot-dashed lines in the top and the middle panels denote 
the results of the self-consistent calculations shown in Fig. 1. 



in Fig. 1 by assuming that the pair potential A(r) is 
proportional to r ■ df /dr, where f{r) is given by Eq. dTTI) . 
We use the proportional constant of —1.107 MeV, that 
leads to the same value for the average pairing gap, 



/Q°°r^dr A(r)p(r) 
r'^dr p{r) 



(17) 



as that in the self-consistent calculations shown in Fig. 
1 for ews(3si/2) = —0.257 MeV. We keep this value in 
varying the depth of the Woods-Saxon potential, — Vq. 

We first keep the particle number to be a constant 
(A^=52) and determine the Fermi energy self-consistently 
within this simplified model. Figure 2 show the results of 
such calculations. For comparison, the top and the mid- 
dle panels also show by the dot-dashed lines the results of 
the self-consistent calculations, which have already been 
shown in Fig. 1. It is remarkable that this model yields 
a similar rms radius to that of the self-consistent calcula- 
tion. The effective pairing gaps show somewhat different 
behavior from those in the self-consistent calculation, es- 
pecially for the pairing gap defined with the canonical 
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FIG. 3: (Color online) The Fermi energy A and the neutron 
number calculated with the simplified HFB model. These 
are plotted as a function of the single particle energy ews 
for the 3si/2 orbit in a Woods-Saxon potential. The solid 
lines are obtained by adjusting the Fermi energy so that the 
neutron number is N—52, while the dashed lines are obtained 
by setting A — ews- 



basis (see the bottom panel). However, it should be em- 
phasized that the pairing gaps stay finite in this calcu- 
lation in the limit of vanishing single-particle energy, as 
in the self-consistent calculation shown in Fig. 1. This 
implies that the self-consistency for the pair potential is 
not important, as far as the rms radius is concerned. 

We next carry out the calculation by setting the Fermi 
energy to be the same as the single-particle energy for the 
3si/2 state, A = e(3si/2)- In this calculation, the number 
of particle changes as we vary the depth of the Woods- 
Saxon potential. The Fermi energy A and the neutron 
number are shown in the upper and the lower panels of 
Fig. 3, respectively, by the dashed lines. For comparison, 
the figure also shows the results of the previous calcula- 
tion shown in Fig. 2, that is, those obtained by adjusting 
the Fermi energy so that the neutron number is a con- 
stant (see the solid lines). The variation of the particle 
number is large, that is, from 47 to 54 in the range of 
single-particle energy shown in Fig. 3. The radii and the 
effective pairing gaps are shown in Fig. 4. As one can 
clearly see, this non-self-consistent calculation yields con- 
siderably different results from the self-consistent calcu- 
lations. Firstly, the reduction of the mean square radius 
of the single-particle orbit is somewhat underestimated, 
although the effect is still large (see the top panel) . Sec- 
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FIG. 4: (Color online) Same as Fig. 2, but by fixing the 
Fermi energy to be the same as the single-particle energy for 
the 3si/2 state, A — ews(3si/2)- 



ondly, the rms radius obtained with this model is com- 
pletely inconsistent with the result of the self-consistent 
calculation as shown in the middle panel. Lastly, the ef- 
fective pairing gaps drops off in the limit of e — >■ (see the 
bottom panel) . Particularly, the lowest quasi-particle en- 
ergy is substantially dismissed, that is a similar behavior 
as that shown in Ref. ll]. Evidently, the claim of Ref. 
[Tlj that the pairing gap disappears in the zero energy 
limit is an artifact of setting A = e. If the Fermi energy 
is determined self-consistently for a given particle num- 
ber, the effective pairing gap persists even if the pairing 
potential is pre-fixed as shown in Fig. 2. 



III. ODD-EVEN STAGGERING OF REACTION 
CROSS SECTIONS 



Let us now investigate how the pairing correlation af- 
fects the reaction cross sections of weakly-bound nuclei. 
To this end, we use the Glauber theory [23l - |28| . In the 
optical limit approximation of the Glauber theory, the 



reaction cross section or can be calculated as |25l-l2S 

2^ 



db 1 



(b) 



(18) 



with 

ix{b) = - j drdr'pp{r)pT{r')TNN{s - s' + b). (19) 

Here, pp and pT are the projectile and the target densi- 
ties, respectively, and b is the impact parameter, s and s' 
are the transverse components of r and r' , respectively, 
that is, s = (r-Bb) ei, and s' = (r'-ef,) et, where = b/b 
is the unit vector parallel to b. Tn]\j is the profile func- 
tion for the A^A^ scattering, for which we assume to take 
a form of fH-il, 



■ NN 



ib) = 



1 



4:71/3 



■ ctnn exp 



213 J 



(20) 



with a^N being the total NN cross section. 

It has been known that the optical limit approximation 
overestimates reaction cross sections for weakly-bound 
nuclei (29l - [33j . In order to cure this problem, Abu- 
Ibrahim and Suzuki have proposed modifying the phase 
shift function x(b) in Eq. ^ to (33} 



ixib) 



1 



drpp{r) 



1 - exp ^- j dr' pT{r')TN]s{s - s' + b) 
i I dr'pTir') 



1 — exp 



drpp{r)rNN{s' - s + b) 



.(21) 



With this prescription, the effects of multiple scattering 
between a projectile nucleon and the target nucleus, and 
that between a target nucleon and the projectile nucleus, 
are included to some extent 33] . 

The upper panel of Fig. 5 shows the reaction 
cross sections for ''4,75,76(^j. _|_ 12q reactions at E=240 
MeV/nucleon, obtained with the phase shift function 
given by Eq. (f2T|) . For the density of the Cr isotopes, 
we use the results of the HFB calculations shown in Fig. 
1. For the density of the target nucleus ^^G, we use the 
same density distribution as that given in Ref. [l^. In 
the actual calculation, we treat the proton-neutron and 
the proton-proton/neutron-neutron scattering separately 
and use the parameters given in Table I in Ref. [33 | for 
the profile function Tnn. In order to evaluate the phase 
shift function, we use the two-dimensional Fourier trans- 
form technique [35j . We give its explicit form in the Ap- 
pendix. The reaction cross sections shown in Fig. 5 show 
a similar behavior as in the rms radii shown in the mid- 
dle panel of Fig. 1, as is expected. That is, the reaction 
cross sections for ''^Gr and ^^Gr are inverted at a small 
biding energy, due to the pairing effect shown in Fig. 1. 
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FIG. 5: (Color online) (The upper panel) Reaction cross sec- 
tions for '■4,75,76qj. _,_ i2(j j-eactions at E=240 MeV/nucIeon, 
obtained with the modified optical limit approximation of the 
Glauber theory, Eq. (I21|l . The density distributions for the 

TA 7^ '7(^ 

' ' Cr are constructed from the HFB calculations shown in 
Fig. 1. (The lower panel) The staggering parameter defined 
byEq. 
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FIG. 6: (Color onhne) The staggering parameter 73 defined 
by Eq. (|22p as a function of the neutron separation energy 
Sn for the odd-mass nuclei. The sohd and the dotted fines 
correspond to the reaction cross sections for 22,23,24q_|_i2q 
and 74,75,76 Qj.^_i2q E=2AQ MeV/nucleon, respectively. 
The dashed and the dot-dashed lines show the results for 
Ne-I-^^C at E=2AQ MeV/nucleon, in which the valence 



30,31,32 



neutron m ' 
tively. 



Ne occupies the 2p3/2 or the If 7/ 2 orbits. 



respec- 



This leads to a large odd-even staggering in reaction cross 
sections for weakly bound nuclei |9|]. 

In order to quantify the OES of reaction cross sections, 
we introduce the staggering parameter defined by 



73 = (-)- 



, aR{A + 1) - 2aR{A) + adA - 1) 



(22) 



where U]i{A) is the reaction cross section of a nucleus 
with mass number A. We can define the same quantity 
also for rms radii. Notice that this staggering parameter 
is similar to the one often used for the OES of bind- 
ing energy, that is, the pairing gap (sgI-Is^. The lower 
panel of Fig. 5 shows the staggering parameter 73 for 
the ''4,75,76Qj. j^mjigj a,s a function of the single-particle 
energy, ews • One can clearly see that the staggering pa- 
rameter 73 increases rapidly for small separation energies, 
and goes up to a large value reaching 73 ~ 80 mb. 

In order to find out a general trend of the stagger- 
ing parameter, Fig. 6 shows the value of 73 for vari- 
ous orbits with 2si/2, 3si/2, 2p3/2 and lf7/2. The values 
for the 2si/2 and 2p3/2 orbits correspond to the reac- 
tion cross sections for the 22,23,24q 30,3i,32]\jg nuclei, 
respectively, calculated in Ref. j^. The value for the 
1^7/2 orbits corresponds to the reaction cross sections for 
30,3i,32]\jg nuclei, obtained with the diffuseness parame- 
ter of the mean-field Woods-Saxon potential of a=0.65 
fm. One can clearly see that 73 for the low I orbits with 
I = and I — 1 show a rapid increase at small separa- 
tion energies, the I = orbits increasing more rapidly 
than the I ~ 1 orbit. In contrast, the high I orbit with 
I = 3 does not show any anomaly in the limit of e — > 0. 
These features are quite similar to the growth of a halo 
structure only in the low I orbits due to a zero or small 
centrifugal barrier. 

The experimental staggering parameters 73 are plot- 
ted in Fig. 7 for Ne isotopes as a function of the neutron 
separation energy for the odd-mass nuclei. We use the ex- 
perimental reaction cross sections given in Ref. Q while 
we evaluate the separation energies with the empirical 
binding energies listed in Ref. [39| . For the neutron sep- 
aration energy for the "^^Ne nucleus, we use the value in 
Ref. [13] ■ The experimental uncertainties of the stagger- 
ing parameter are obtained as 



(573 



^[SaniA + 1))2 + A{5aR{A)f + {5<Jn{A - l)f 



(23) 

where 6afi{A) is the experimental uncertainty for the re- 
action cross section of a nucleus with mass number A. 
The figure also shows by the dashed line the calculated 
staggering parameter for the 30,3i,32]\jg nuclei with the 
2p3/2 orbit, that has been shown also in Fig. 6. One sees 
that the experimental staggering parameter agrees with 
the calculated value for ^0,31,32^,^^ nuclei when one as- 
sumes that the valence neutron in ^^Ne occupies the 2p3/2 
orbit. Furthermore, although the structure of lighter 
odd-A Ne isotopes is not known well, it is interesting to 
see that the empirical staggering parameters closely fol- 
low the calculated values for the 2p3/2 orbit. This may 
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FIG. 7: (Color online) The experimental staggering parame- 
ter 73 of reaction cross sections defined by Eq. ((22} for the Ne 
isotopes with the ^^C target at £^=240 MeV/nucleon. This 
is plotted as a function of the neutron separation energy S„ 
of the odd-A nuclei. The experimental data for the reaction 
cross sections are taken from Ref. Q, while the empirical 
separation energies are taken from Refs. [3^.|40|]. The dashed 
line is the calculated staggering parameter for the 30,3i,32-|\^g 
isotopes, assuming that the valence neutron in of ^^Ne occu- 
pies the 2p3/2 orbit. 



indicate that the low-Z single-particle orbits are appre- 
ciably mixed in these Ne isotopes due to the deformation 

effects mm- 



IV. SUMMARY 

We have studied the odd-even staggering (OES) of 
the reaction cross sections by using the Hartree-Fock- 
Bogoliubov model. To this end, we have introduced the 
staggering parameter 73 defined with a 3 point difference 
formula in order to clarify the relation between the mag- 
nitude of OES and the neutron separation energy. We 
have shown that the OES parameter increases largely for 
low-/ orbits with I = and / = 1 at small separation 
energies. The experimental staggering parameter for the 
Ne isotopes show a similar increase. On the other hand, 
we have found that the staggering parameter stays al- 
most a constant value, 73 '~ 2 mb for higher I orbits with 
e.g., I = 3. The increase of 73 is induced by the finite 
pairing correlations in the zero separation energy limit. 
In this respect, we have shown that the effective pair- 
ing gap for the 3si/2 orbit in the ^^Cr nucleus persists 
even in the limit of vanishing separation energy. This 
remains the same even if the pair potential is prefixed, 
as long as the chemical potential is adjusted to keep the 
particle number to be the same. We have shown that 
such simplified HFB model well reproduces the results of 
the self-consistent HFB model for the root-mean-square 
radius. 

The staggering parameter proposed in this paper pro- 
vides a good measure for the OES of reaction cross sec- 
tions. Further systematic experimental studies would 



be helpful in order to clarify the pairing correlations in 
weakly-bound nuclei and in the limit of zero neutron sep- 
aration energy. 
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Appendix A: Evaluation of phase shift function with 
the Fourier transform method 



In this paper, we evaluate the phase shift functions 
given by Eqs. (jl9p and (1211) using the Fourier transform 
technique [35]. First we notice that Eq. (fT9|) can be 
expressed as 

ix{b)^- j dsds'p'^;\s)p^^\s')TNN{s-s' + b), (Al) 
with 

Pp\s)= Jdzpp{s,z), pi^\s') = J dzpT{s',z). 

The two-dimensional Fourier transform of ix{t>) in Eq. 
(IT9t then reads 



tx{q) = J dbix{b)e^'i-^, (A3) 
= - J dbe'K^+^-^'^TNNib + s-s') 
X J dse-'1-%f{s) I ds' p'>^\s'){M) 

where T]^pf,pp, and are the two-dimensional 

Fourier transform of TNPf,p^p\ and p^\ respectively. 
For the profile function given by Eq. ((20)) . its Fourier 
transform reads. 



NN[q) = , „ CTNN ■ 2p Trexp 
47rp 



(A6) 



The Fourier transform of the density distribution is eval- 
uated as 

p(^)(q) = J dzdse'^i-^ pis,z)^piQ), (A7) 

where p{Q) is the three-dimensional Fourier transform of 
the density at Q = {q,Qz = 0). For a spherical density. 
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p(r), p{Q) depends only on |Q| = q, that is, 

/•OO 

p'-'Hq) = Kq) - 4^ / r^drp{r)jo{qr), (A8) 



where joiqr) is the spherical Bessel function of zero-th 
order. Taking the inverse Fourier transform of Eq. (jASp . 
the phase shift function is calculated as 



ixib) 



dq 

(2^ 



fjviv(g)pp((z)pT(<z)e-'9-b,(A9) 



§ VMN{q)~Pp{q)PT{q)Mqh)iklQ) 



where Jo{qb) is the Bessel function of zero-th order. A 
similar technique has been used to evaluate a double fold- 
ing potential in heavy-ion reactions [43l - l45| . 

One can apply the same method to evaluate the phase 
shift function given by Eq. (PT|) . First notice that 

TNT{b) = l-exp(^- J dr'pT{r')rNN{b-s')yAn) 

= l-exp^-^ ^tNN{q)pT{q)Jo{qb)^ , 

(A12) 



depends only on 6 = \b\. This leads to 



ix(b) 



in 



qdqfNT{q)pp{q)Joiqb) 



— / qdqfNp{q)pT{q)Mqb), (A13) 



where f ArT('?) is calculated as 



rjvT(<z) = 2^/ bdbrMTib)Mqb). (A14) 

JO 
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